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Abstract 

Using 3.3 pb" 1 of data collected with the CMD-2 detector in the 720 - 840 MeV 
cm. energy range, the branching fraction of the conversion decay uj — > 7r°e + e~ has 
been measured: B(lo — ► 7r°e + e~) = (8.19 ± 0.71 ± 0.62) • 10~ 4 . The upper limits for 
the branching fractions of the following conversion decays have been obtained at 
the 90% confidence level: B(p -> ir°e + e-) < 1.6 • 10~ 5 , B(p -> r/e+e") < 0.7 • 10" 5 
and B(u — * -ne + e~) < 1.1 • 10~ 5 . 
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PACS: 13.40.G, 13.40.H, 14.40.E 



1 Introduction 

Measurement of branching fractions and transition form factors of conversion 
decays provides an important test of vector dominance model [1,2] and an 
accurate background estimation in searches for quark-gluon plasma involving 
a lepton pair [3,4]. 

The expected branching fractions for the p and u conversion decays are B(u> — > 
7r°e + e-) = (7.2 - 8.0) • 10" 4 , B(p -> 7r°e+ e -) = (4.1 - 6.5) • 1(T 6 and an 
order of magnitude smaller for decays with muons and/or 77 meson in the final 
state [5,6,7]. Various conversion decays of the <fi meson were studied at CMD- 

2 [8,9] and SND [10,11]. For the p and u mesons experimental information is 
rather scarce: only the branching fractions for the decays uj — > 7r°e + e~ and 
uj — > n°p + p~ were measured in experiments at ND [12] and Lepton-G [13] 
with an accuracy ~ 30%. 

Conversion decays of a vector meson V into a pseudoscalar meson P and lepton 
pair l + l~ (V — > Pl + l) are closely related to the corresponding radiative 
meson decays V — > P7 [5] and provide a possibility to study a transition 
V — > P form factor, Fvp(q 2 ), as a function of squared mass of virtual photon, 

q 2 = MUin-). 

In this work we present the result of the studies of the p and uj conversion 
decays into a e + e~ pair and pseudoscalar meson (n or rj) performed with the 
CMD-2 detector at the VEPP-2M collider [14]. The analysis is based on a 
data sample collected at 19 energy points in the 720 - 840 MeV cm. energy 
range and corresponding to 3.3 pb _1 of integrated luminosity. This statistics 
contains about 3.3-10 6 p and 1.8-10 6 u decays. More detailed description of 
this analysis can be found in [15]. 



2 Experiment 

The general purpose detector CMD-2 has been described in detail elsewhere [16] 
The tracking system consists of the cylindrical drift chamber (DC) with 250 p 
resolution transverse to the beam plane and double-layer multiwire propor- 
tional Z-chamber, both also used for the trigger. The tracking system is placed 
inside a thin (0.38 X ) superconducting solenoid with a field of 1 T. The bar- 
rel Csl calorimeter with a thickness of 8.1 X placed outside the solenoid has 
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energy resolution for photons of about 9% in the energy range from 100 to 
700 MeV. The angular resolution is of the order of 0.02 radians. The end-cap 
BGO calorimeter with a thickness of 13.4 X placed inside the solenoid has 
energy and angular resolution varying from 9% to 4% and from 0.03 to 0.02 
radians, respectively, for the photon energy in the range 100 to 700 MeV. The 
barrel and end-cap calorimeter systems cover a solid angle of 0.92 x An radians. 



3 Data analysis 

The decay uj — > 7r°e + e~ has been studied using the n° dominant decay mode 
7T° — > 77. It corresponds to a final state with two opposite charge particles 
and two photons. 

One of the significant resonant backgrounds comes from the uj — > 7r + 7r~7r° de- 
cay which has the same topology of the final state and more than three orders 
of magnitude larger probability. Another source of resonant background is the 
uj — > 7r°7 decay followed by the Dalitz decay of the 7r° or 7-quantum conversion 
in the material in front of the drift chamber. Most conversions occur in the 
beam pipe, at the distance of 1.8 cm from the interaction point. Since the DC 
spatial resolution is not sufficient to separate events with conversions in the 
beam pipe from those where it occurs at the interaction point, the contribution 
of this background was subtracted based on Monte Carlo simulation (MC). 
The non-resonant background includes contributions from the following QED 
processes with the same final state topology: e + e~ — > e + e~77, e + e~ — > 37 fol- 
lowed by 7-quantum conversions, e + e~ — > e + e~7 with one background photon 
as well as two-quantum annihilation followed by a 7-quantum conversion and 
one background photon in calorimeters. 

The data analysis consists of two parts: a measurement of the branching frac- 
tion and a study of the transition form factor. The selection criteria used for 
the measurement of the branching fractions were tuned to select a "pure" set 
of events under study and thus reject events with a large invariant mass of 
e + e~ which contain a lot of background. However, this very range is of main 
interest for a study of the transition form factor. Therefore, for the latter a 
special set of selection criteria including a technique of e/n separation [9,17] 
was applied to suppress the main background from 3ti events. 

3.1 Selection of p{uj) — > 7r°e + e~ events 

At the first stage of selection, the following criteria were applied to events 
with two tracks and at least two photons to enrich a data sample with events 
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of the studied decays: 



• the photon energy E 1 > 40 MeV and its polar angle 0.5 < 7 < ir — 0.5 to 
suppress background photons in the calorimeters; 

• the impact parameter of the tracks p < 1 cm and Z-coordinate of the vertex 
\Z ver t\ < 5 cm to reject cosmic rays and beam background events; 

• the total momentum of the tracks p = \P\ + P2I does not strongly differ 
from the photon momentum p 7 in the U) —>■ 7r°7 decay at a given energy 
\P ~ P-y\ < 35 MeV/c to suppress uj — > 7r + 7r~7r° events as well as uj — > 7r°7 
events followed by the Dalitz decay of the tt°. This condition is illustrated 
with horizontal lines in Fig. 1 where events satisfying previous requirements 
are shown; 

• the angle between the total momentum of the tracks and each photon is 
more than 1.7 to suppress QED events; 

• the invariant mass of the electron-positron pair and most energetic photon 
M inv (e + e~^i) is less than 1.9 • E beam to suppress e + e~ — > 77 events followed 
by a conversion of one 7-quantum; 

• the opening angle of the tracks Aip < 0.5 to suppress events of the uj — ► 
7r + 7r - 7T° decay, this condition is shown in Fig. 1 with a vertical line. 

Using these criteria 390 events were selected. 




Fig. 1. Opening angle of two tracks Aip 
versus their total momentum p. The 
horizontal lines show the cut on the 
total momentum and the vertical one 
shows the cut on the opening angle of 
the tracks. 



Fig. 2. The Miny^jj) distribution (his- 
togram). The results of the fit are shown 
with a solid (signal) and dotted (back- 
ground) line. 



At the next stage, a fit of the 77 invariant mass distribution {M inv ( / y / y)) was 
performed at each energy point to determine the number of events with tt° in 
the final state denoted hereafter as "7r°2c" (uj — ► ir°e + e~ , uj — ► 7r°7 followed by 
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a 7-conversion and uj — > 7r + 7r~7r°). The parameterization included a logarith- 
mic Gaussian for the signal and a second order polynomial for the background. 
The shape of these functions was fixed, so the floating parameters were the 
numbers of events of each type only. The shape of background was taken from 
a study of the M inv {^^) distributions from the uj decays into 7r°7 and n + n^n° 
specially selected in the whole energy range. The shape of the signal was ob- 
tained from the fit of the Mj m) (77) distribution over all selected events with 
floating parameters of the function describing the signal, see Fig. 2. The total 
number of selected 7r°2c events is 316. 



3.2 Study of main background processes 

At this stage of analysis the expected number of background events was eval- 
uated using the Aip distribution. Its difference for events of the process under 
study and those of background was employed in the analysis: the former are 
strongly peaked near zero while events coming from the uj — > 7r + 7r~7r° decay 
have a smooth wide distribution. Events of the uj — > 7r°7 decay with conver- 
sion, which are also strongly peaked near zero, predominantly populate the 
angular range Aip < 0.3. The Aip distribution for the data after imposing all 
selection conditions but that on Aip is shown by crosses in Fig. 3. Using the 
modified requirement 0.3 < A^ < 0.9 we reject most of the events coming 
from the uj — > 7T°7 decay with conversion and keep events under study as well 
as background events from the uj — > 7r + 7r~7T° and QED. The MC Aip distribu- 
tion for the main background of uj — > 7r + 7r~7r° shown with a hatched histogram 
in Fig. 3 was approximated with a function f(x) ~ x a with a = 1.63 ±0.10 in 
the Aip < 0.9 angular range. The energy dependence of the number of events 
selected at the first stage (with modified requirements on Aip) was used to 
determine the number of uj — > 7r + 7r~7r° events. 

To this end, the number of selected events at each point was written as 

Ni = (cr res (si) (1 + Si) + o-QED(si)) ■ U, (1) 

where ctqed(s) = °~qed ' m t/ s -> °Ves(s) is the cross section of the uj meson 
production at the squared cm. energy s, Li denotes the integrated luminosity 
at the i-th energy point, and 5i is a radiative correction. The floating param- 
eters are o-q ED and of es — cross section values at the resonance peak. The 
resulting number of resonance events was calculated as a sum over all en- 
ergy points: iV res = I]o" r cs(s«) Li (1 + Si) = 210 ± 17. After subtraction of the 
uj — > 7r°e + e~ contribution, the number of 7r + 7r~7r° events N n + n - n o = 180 ± 20 
in the 0.3 < Aip < 0.9 range. The above description of the Aip distribution 
was used to calculate the number of uj — > -n + n~n° events in the Aip < 0.5 
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Fig. 3. The Aip distribution for the data (crosses) and u — ► 7T + 7T~tt simulation 
events (hatched histogram). All selection criteria but that on Aip were applied. The 
number o/tt+tt^tt MC events was approximated with a function x a shown with a 
solid line. 

range: 



N 37r = 40.6 ± 6.4 ± 4.6. 



(2) 



The systematic error comes from the uncertainty of the Aip distribution ap- 
proximations and that of the subtraction of the ui — > 7t°e + e~ contribution. 
Since events of u — > tt°7 decays could not be distinguished from those under 
study, the detection efficiency for the former was determined as e^ +e ■ P C onv, 
where e^ +e denotes the detection efficiency for 7T°e + e~ events and P con v is 
a probability for a monoenergetic photon from the u) —>■ tt 7 decay to convert 
in the material in front of DC. The P C onv value was obtained from simulation: 



Pr, 



1.91 ±0.06) • 10" 



(3) 



The error of the P con v value above comes from the uncertainties in a thickness 
and composition of materials in front of DC. 



3.3 Approximation o/tt 2c events 



At the last stage of analysis, the energy dependence of the number of tt 2c 
events is fitted with a sum of the contributions from the u — > 7T°e + e~ , tt + 7t~7t°, 
tt°7 decays and possibly remaining QED events: 



-^r°2c,i = N n e+e - 



+ N^ i + N 37r>i + N Q ED, 



(4) 
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The number of n°e + e events at the i-th energy point was described by the 
expression (5): 

A^e+e-.i = a 7r o e+e -(s i ) (1 + 5i)£dcM e e 'trig, it3(n° ->77), (5) 

where Lj is an integrated luminosity, <5j is a radiative correction, e^lt^ , £Ai/>,i, 
£tri g ,i denote the detection efficiency, efficiency of the reconstruction of close 
tracks and trigger efficiency at the i-th energy point, respectively. The en- 
ergy dependence of the Born cross section a n o e + e - (s) was written using the 
relativistic Breit-Wigner approach with the p and u meson contributions: 

Vn°e+e-(s) = ' \ A P ( S ) + A M) + a of • (6) 



Here q(s) is a phase space factor: 

?(«) = ~y I 1 - -j-J> ( 7 ) 



-Ay(s) is an amplitude of the vector meson V: 



v 



D v (s) q^(m v Y 



(8) 



and the additional constant ao in the amplitude describes a possible contri- 
bution of higher resonances. The quantities my and a v are the vector meson 
mass and the cross section at its peak, respectively, with a v calculated without 
taking into account other contributions, and fy = e l ^ v is a phase factor. To 
describe the p-ou interference, the model with energy-independent interference 
phases was used with W set to zero and 4> p = —13° [18]. The l/D v (s) is a 
vector meson propagator described by the expression: 

D v (s) = m v — s — im v T v (s). (9) 



The quantities Ty(s) and IV = Ty{m v ) are vector meson width at the squared 
cm. energy s and at the vector meson mass, respectively. 

The number of 7r°7 events was written as 

N*°y,i = °"7r° 7 ( S U (! + Si 7 ) £ dlCi C Pconv £ AV.,^trig,i B{ir° -> 77), (10) 



with <7„-o 7 (s) taken from [18] and a photon conversion probability P con v taken 
from (3). To describe the energy dependence of ir + iT~ir events, the relativistic 
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Breit-Wigner with the u meson contribution was used: 

N 3 n,i = VBw{Si) U (1 + bf) e 3n . (11) 



Since the total number of 7r + 7r~7r° events was previously found, see expres- 
sion (2), the normalization factor e 37r was determined from the sum N 3n = 
S c"Bw( s i) Li (1 + 5f n ) e 3n . A possible contribution of QED events which can, 
e.g., appear because of the incorrect background subtraction was evaluated 
from the following expression: 

in? 

Nqed, i = ctq ED U — . (12) 



The detection efficiencies e^l^ e were determined using the Monte- Carlo sim- 
ulation taking into account initial and final state radiation. The calculation of 
radiative corrections 5 followed Ref. [19]. 

The events under study can be triggered by two independent triggers: charged 
and neutral ones. The overall trigger efficiency is estimated to be 99%. The 
estimation of this value has been performed by analysis of the trigger signals 
in selected events. 



Since Monte-Carlo simulation does not completely describe the experiment, 
a correction e^ip for a difference between the efficiencies of close track recon- 
struction in simulation and experiment was included to describe the data. 
Its value was obtained using events of uj — > 7r + 7r~7r° decays followed by the 
conversion decay n° — > e + e~7 with a similar distribution. The integrated 
luminosity was determined using events of large angle Bhabha scattering with 
radiative corrections taken into account according to [20]. 

Table 1 shows the detailed information about experiment including efficiencies 
of detection, trigger and close track reconstruction together with the number of 
7r°2c events and visible cross section calculated from the following expression: 

/V cxp 

vis _ 1 y Ti-o 2c (1^ 

^ ~L(1 + S) £5e°t e+e " ^ ^trig B{lfi - 77 ) ■ 1 ' 



The energy dependence of the number of 7r°2c events was fitted using the 
maximum likelihood method. The minimization function is shown here: 

c = ^ ( ( N ^2 C ,i-K7 2Cti )\ (iy 37r - jv 37r ) 2 K- m j 2 (f „ - r j 2 
t-S °h } n** e(r.) 2 ; 
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Table 1 

The energy, integrated luminosity, radiative correction, detection efficiency, the re- 
construction efficiency of close tracks, trigger efficiency, number of 7r°2c events and 
visible cross section of 7r°2c events. 



v/s, MeV 


L, nb- 1 


5 


£ det 




£trig 




rr vis nh 
°>2c> 110 


720 


183.86 


-0.099 


0, 


149 


0.906 


0.989 


o-8±5:| 


04 +0 - 06 


750 


148.66 


-0.119 


0. 


149 


0.906 


0.989 


1 5 +L6 
1 -°-0.9 


09 +0 - 10 

u - uy -0.06 


760 


157.11 


-0.152 


0. 


.149 


0.906 


0.989 


q 4+2.2 
°-^-1.5 


20 +0 - 13 

u - zu -0.09 


770 


45.81 


-0.186 


0. 


.149 


0.906 


0.989 


o.ol5:8 


o.oolg:£ 


774 


119.57 


-0.211 


0. 


.148 


0.906 


0.989 


3-Hfi 


0.26t°;g 


778 


128.25 


-0.218 


0. 


.148 


0.906 


0.989 


19.433 


1 ™+ - 39 


780 


132.82 


-0.229 


0, 


.148 


0.906 


0.989 




1 61 +0 - 38 


781 


187.09 


-0.224 


0. 


148 


0.906 


0.989 


28.8+jJ 


i c; 9 +0.31 


782 


240.98 


-0.211 


0. 


.148 


0.906 


0.989 


56.8l 8 7 | 


2 30+ 34 
Z - OU -0.32 


783 


187.34 


-0.191 


0. 


.148 


0.906 


0.989 


54.i±f:§ 


9 77 +0.42 
z -' '-0.39 


784 


237.05 


-0.167 


0. 


149 


0.906 


0.995 


56.3±fi 


9 on+0.33 
z - zu -0.31 


785 


191.51 


-0.138 


0. 


149 


0.906 


0.998 




1 63+ 035 
i.oo_ 26 


786 


108.63 


-0.112 


0, 


149 


0.906 


0.998 


9.2t 3 2 A 7 


71 +0 - 26 

u - ' i -0.21 


790 


117.03 


0.053 


0. 


149 


0.906 


0.998 




°- 60 -o:il 


794 


121.81 


0.260 


0. 


.149 


0.906 


0.998 


5-i±?:g 


27+ - 14 
u - z/ -0.10 


800 


195.50 


0.427 


0. 


.150 


0.907 


0.998 




o n7 +0 - 06 

u - u '-0.04 


810 


181.17 


0.828 


0. 


151 


0.907 


0.998 


2 6 +2 -° 
z -°-1.3 


o n7+ - 05 

u - u '-0.03 


820 


185.52 


1.023 


0. 


151 


0.907 


0.998 


o 7+2.0 
z -'-1.3 


06+ - 05 

u - UD -0.03 


840 


457.87 


1.260 


0. 


151 


0.907 


0.998 


O + 2.1 

°- z -1.5 


03+ - 02 



where i\^o 2c ,i = N n o e+e - :i + iV„.o 7ji + N 377ji + N QED:i and iV^ Cji is the experi- 
mental number of 7r°2c events at the i-th energy point. The quantities 4,j are 
asymmetric variances of the number of 7r°2c events. The floating parameters 
are two branching fractions B(uj — > 7r°e + e~) and B(p — > 7r°e + e~), and the con- 
stants ao and Oq^. Since the statistics is not sufficient to determine precisely 
main parameters of the ui meson, the fit allows uj meson mass m w , width 
and the number of 7r + 7r~7r° events N 3n (three last terms in expression (14)) to 
float around their expected values within their estimated uncertainties. The 
central value and error and were taken from PDG [21] while those for 
N 3n were taken from (2). The p meson parameters were fixed from the work 
[22]. 
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Table 2 

The fit results in various models. Only statistical errors are shown. 

Model B(u), 10" 4 B(p), 10" 4 a QED , pb a , nb 1 / 2 x 2 /n.d.f. 

I 8.19±0.81 fixed -0.1±2.5 = 19.11/17 

II 8.20±0.77 fixed = 0.10±0.14 18.17/17 

III 8.19±0.71 fixed =0 =0 19.11/18 

IV 8.33±1.37 0.032±0.081 = = 19.06/17 

Four different models were used. The ratio of o"° and a Q p values, cross sections 
at the p and uj meson peak, respectively, was fixed in some models from the 
results of the study of the process e + e~ — > [18], so that 

a°(p -> 7r°e+e-) _ a°(p -> 7r° 7 ) nA n _ 3 



(T°(c<j — > 7r°e+e ) cr (o; — > 7r°7) 



= (3.7 ± 1.0) • 1(T'\ (15) 



The theoretical precision of the approximate equality in (15) is evaluated as 
3% for the applied selection criteria because of small q 2 in selected events. 
The expression (15) fixes the ratio between the branching fractions B(p — > 
ir°e + e~)/B(uj -> 7r e + e -) = (1.6 ± 0.4) • 10" 3 . 

A fit in the first model was performed with floating <Jq ED to evaluate a possible 
contribution of QED events. As a result, the contribution of QED to 7r°2c 
events is negligible, so in the other models this parameter was fixed to zero. 
Later, the fit in model II had a floating a to evaluate a possible contribution 
from higher resonances. As a result, this contribution a does not differ from 
zero within a statistical error. In model III a fit with a fixed ratio cr°/<r° was 
performed to obtain a branching fraction for the uj — > 7r°e + e~ decay 

B(uj -> n°e + e-) = (8.19 ± 0.71 ± 0.62) • 10~ 4 . (16) 



The graphical demonstration of this fit together with the visible cross section 
of the 7r°2c events is shown in Figure 4. As a result of this fit the number of 
uj — > 7r°e + e~ events was determined to be 232. Finally, a fit in model IV with 
floating B(uj — > 7r°e + e~) and B(p — > 7r°e + e~) was performed to set an upper 
limit on the branching fraction for the p — > n°e + e~ decay: 

B(p -> 7i°e + e-) < 1.6 • 10~ 5 (90% C.L.). (17) 



The obtained value of the uj — > n°e + e~ branching fraction (expression (16)) 
agrees with the previous measurement from ND [12] but has two times better 
accuracy. It is also consistent with the theoretical predictions [6,7]. 
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Cross-section fit by Breit-Wigner 




10 b-l— L_L_X L 4 — , — ;— y- — f — ,— y- -;- — j — ,— y- T J - C I "f H " 

720 740 760 780 800 820 840 

is, MeV 

Fig. 4. Visible cross section of 7r° '2c events. The fit in model III is shown. 
3.4 Systematic errors 

The main sources of systematic uncertainties of the branching fractions are 
listed in Table 3. A contribution of each parameter to the total systematic er- 



Table 3 




The main sources of systematic uncertainties 




Source 


Uncertainty, % 


Reconstruction efficiency of close tracks £a</> 


4.7 


Background subtraction 


3.6 


Detection efficiency £^ +e 


2.2 


Trigger efficiency e t rig 


2.0 


Parameters of p and to 


2.0 


Mixing parameter a = a^/a® 


1.9 


Integrated luminosity 


1.4 


Form factor model 


1.2 


Radiative corrections 


1.2 


Probability of conversion 


1.0 


Total 


7.6 



ror was estimated as a change of the branching fraction when this parameter 
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was varied within its measurement uncertainty. The main contribution comes 
from the reconstruction efficiency of close tracks and it is completely deter- 
mined by statistics of the test samples of the uj — > 7r + 7r~7r° decays followed by 
7T° — > e + e~7. The background subtraction error originates from the inexact 
knowledge of the shape of the M{^) distribution for background events and 
from the error in the number of 7r + 7r~7r° events. Since the reconstruction of 
close tracks is sensitive to the presence of the charged trigger, a contribution of 
the trigger efficiency was conservatively evaluated as a change of the branching 
fraction when we additionally demanded that the charged trigger was on. The 
detection efficiency error was determined by varying selection criteria. The 
uncertainty due to the determination of integrated luminosity comes from the 
selection criteria of Bhabha events, radiative corrections and calibration of the 
DC and Csl calorimeter. The uncertainty of the model of the transition form 
factor was evaluated by a difference in the detection efficiency when simple 
VDM with a single p meson as a transition particle and generalized VDM 
(p and p') were used in simulation. The uncertainty of radiative corrections 
comes from the dependence on the emitted photon energy and the accuracy 
of theoretical formulae. The resulting systematic uncertainty of the branching 
fraction quoted in Table 3 is 7.6%. 



3.5 Search for r]e + e events 

A search was performed in three main decays modes of the 77 meson: 3n°, 
7T + 7r~7r and 27. Selection criteria for a e + e~-pair were the same as for the 
7i°e + e~ final state: p < 1 cm and \Z vert \ < 5 cm; the difference between the 
total momentum of the tracks p = \p± +P2I and the momentum p 7 in the 
uj — > 777 process at a given energy \p — p 7 | < 30 MeV/c; the opening angle 
of the tracks Aip < 0.5. In the 77 — > it + it~it decay mode with a four-track 
final state a e + e~-pair was identified as a pair with opposite charges and the 
smallest opening angle. The following cuts on photons were applied: 

• 77 — > 3tt°: the photon energy threshold was decreased to 30 MeV because of 
a higher multiplicity of soft photons in the final state, N 7 > 5; 

• 77 -> tt+tt-tt : N 1 = 2, M inv (e+ , e" , 71) , M inv (e+ , e~ , 72) > 160 MeV/c 2 to 
suppress events from decays uj — > 7r + 7r~7r° followed by n° — > e + e~7; 

• 77 — > 27: N y = 2, the soft photon energy i? 7)2 > 175 MeV to suppress QED 
events, M inv (e + , e - , 71), M inv (e + , e~, 72) > 200 MeV/c 2 to suppress events 
from decays uj — > 7r°7 followed by n° — > e + e~7, the total momentum of the 
tracks and photons P = \pi + p 2 + p 7 ,i + p 7 ,2| < 150 MeV/c. 

The main background for the i] — > 7r + 7r~7r° and 77 — > 27 decay modes comes 
from the QED processes while for the 77 — > 37r° decay mode it comes from 
the decay uj — > 7r°7r°7 with a "fake" photon in calorimeters. The result of the 
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Table 4 

The branching fraction, the detection efficiency at s = m^, the reconstruction effi- 
ciency of close tracks, the trigger efficiency, the number of selected events and the 
average number of background events for analyzed rj meson decay modes. 



Decay mode 


B, % 


e det 




^trig 


-Nexp iV back 


T) -^3vr° 


32.51±0.29 


0.070 


0.89 


0.98 


< 0.1 


Tj — > 7T + 7T~7r 


22.6±0.4 


0.018 


0.89 


0.99 


0.2 




39.43±0.26 


0.059 


0.89 


0.98 


3 4.9 



analysis is shown in Table 4. 

The total number of selected events for each 77 decay mode was expressed as: 

n 

N = ^2<TBw(si) Li(l + bf 6 )e det (si)eA^e U igB(r] ^ final). (18) 



The cross section parameterization is similar to formula (6) and includes rel- 
ativists Breit-Wigner contributions of the p and u mesons, but without the 
do term. The interference model with a constant phase §^-p = is used. We 
used the results of our study of the process e + e~ — > 777 [23] to fix the val- 
ues of the amplitude: a°(p — > r]e + e^)/a°(uj — > r\e^e ) = 0.43 ± 0.09, from 
which the following relation between the branching fractions can be obtained: 
B(p — > rje + e~)/B{uj — > i]e + e~) = 0.65 ± 0.14. Since the detection efficiency 
strongly depends on s, its value was calculated at each energy point. The ra- 
diative corrections S ve+e , trigger efficiency e t ri g as well as the reconstruction 
efficiency of close tracks were calculated by the same method as for the 
7r°e + e~ events. The average number of background events has been evaluated 
by MC simulation. The simulation of corresponding events included generation 
of "fake" photons [24] . The results of all analyzed 77 decay modes are taken into 
account to calculate the upper limits on the branching fractions of the p and 
to meson decays. The total number of selected events has been described by 
summation over expressions (18) applied to each r\ decay mode. To calculate 
upper limits, the Feldman-Cousins approach [25] was applied with 3 events 
observed and 5.1 background events expected. The final values of the upper 
limits were increased by 13% in accordance with the systematic error. The 
obtained upper limits exceed theoretical predictions by a factor of 2 [5,6,7]: 



B{p -> 7]e + e-) < 0.7 • 10~ 5 (90% C.L.), (19) 
B{uj -> r]e + e-) < 1.1 • 10" 5 (90% C.L.). (20) 
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3.6 Study of the ujtt transition form factor 



The electromagnetic transition form factor F(q 2 ) provides information on the 
electromagnetic structure of interacting particles. Its dependence on squared 
invariant mass of virtual photon q 2 can be approximated in a small q 2 range 
as 

F(q 2 ) = l + b.q 2 , (21) 



where the parameter b is a slope of the transition form factor. In experiment 
the information about photon virtuality in events of uj — > 7T°e + e~ decay could 
be taken from the invariant mass of the lepton pair Mi nv (e + e~) = q 2 . Since 
the previously used selection criteria leave events with small q 2 only and the 
detection efficiency strongly depends on M inv (e + e~), the cuts on total momen- 
tum, angle between the total momentum of tracks and photons, M im) (e + e~7i) 
were omitted and the cut on was increased up to 2.5 radians. To suppress 
background mainly originating from uj — > 7r + 7r~7r°, the kinematic fit requiring 
energy-momentum conservation and the procedure of e/ir separation based on 
the energy deposition of charged particles in calorimeters [9,17] were addition- 
ally used. 

The selected events were divided into groups of M inv (e + e~). The number of 
uj — > 7r°e + e" events in each group was determined from the analysis of the 
77 invariant mass distribution. After that the q 2 dependence of the obtained 
number of events was fitted to the following function from [5] 



dN n gl { 4mh 1/2 /, 2m f 2 > ± 



x 



q 2 \ 2 4m 2 g 2 

1 + 



ml -ml J (m 2 -m 2 ) 5 



3/2 

,2\|2 



\F^(q 



where a is the fine structure constant. The detector resolution on invariant 
mass a q = 15 MeV/c 2 was taken into account. The floating parameters are 
the normalization factor A and the slope of the transition form factor b. The 
result of the fit is shown in Fig. 5, the slope value b = 2.5 ±3.1 GeV~ 2 . 
The total systematic error is estimated to be 10%. The main contributions 
to it are the e/ir separation procedure and selection criteria. This error is 
smaller than a statistical one for all invariant mass groups except for the group 
with q < 50 MeV/c 2 . The obtained value agrees with the VDM prediction 
b = l/m 2 = 1.7 GeV~ 2 within statistical errors. 

It is worth noting that modes with fj, + fi~ in the final state are more convenient 
for studies of the transition form factor \F VP (q 2 )\ [13]. Unfortunately, their 
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q, invariant mass of e e -pair, GeV 

Fig. 5. Results of the fit of the data (dots with error bars) on the electromagnetic utt 
transition form factor \F(q 2 )\ 2 (solid line). The VDM form factor is shown with a 
dotted line. 

study at CMD-2 is complicated because of the large background coming from 
the dominant decay mode of the to meson: uj — > 7r + 7r _ 7r° and insufficiently 
powerful h/tc separation at low momenta. 



4 Conclusions 

Using a data sample of 3.3 pb _1 in the 720 - 840 MeV cm. energy range, the 
branching fraction for the u — > 7c°e + e~ decay has been obtained: 

B(uj -> n°e + e-) = (8.19 ± 0.71 ± 0.62) • 10" 4 

and the following 90% C.L. upper limits were set: 

B{p -> n°e + e-) < 1.6 • 10~ 5 , 
B(p^T]e + e-)<0.7- 10" 5 , 
B(u -> r]e + e~)< 1.1- 10" 5 . 
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